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ABSTRACT 

We present first results from a high resolution multi-band survey of the Westerlund 2 region with 
the Hubble Space Telescope. Specifically, we imaged Westerlund 2 with the Advanced Camera for 
Surveys through the F555IU, F814IU, and F658N filters and with the Wide Field Camera 3 in the 
F125IU, E160IU, and F128N hlters. We derive the first high resolution pixel-to-pixel map of the 
color excess E{B — V)g of the gas associated with the cluster, combining the Ha {F658N) and Pa/3 
(^128^) line observations. We demonstrate that, as expected, the region is affected by significant 
differential reddening with a median of E{B — V)g = 1.87 mag. After separating the populations 
of cluster members and foreground contaminants using a (F814IU — E160IU) vs. F814IU color- 
magnitude diagram, we identify a pronounced pre-main-sequence population in Westerlund 2 showing 
a distinct turn-on. After dereddening each star of Westerlund 2 individually in the color-magnitude 
diagram we find via over-plotting PARSEC isochrones that the distance is in good agreement with 
the literature value of 4.16 ± 0.33 kpc. With zero-age-main-sequence fitting to two-color-diagrams, 
we derive a value of total to selective extinction of Rv = 3.95 ± 0.135. A spatial density map of the 
stellar content reveals that the cluster might be composed of two clumps. We estimate the same age 
of 0.5-2.0 Myr for both clumps. While the two clumps appear to be coeval, the northern clump shows 
a ~ 20% lower stellar surface density. 

Subject headings: techniques: photometric - stars: early type - stars: pre-main sequence - HII regions 
- open clusters and associations: individual (Westerlund 2) - infrared: stars 


1. INTRODUCTION 

Wi th a stellar mass of M > 10^ Mq (jAscenso et al.l 
1200711 and an estimated molecular cloud mass o f (1.7 ± 
0.8 - 7.5) X 10® Mq (iFurukawa et al.l l2009t iDamd 
120071 both based on millimeter CO spectroscopy), the 
young, massive sta r cluster Westerlund 2 (hereafter Wd2; 
IWesterhmdl 1196111 is one of the most massive young 
star clusters known in the Milky Way (MW). It is 
embe dded in the H II region RCW 49 (iRodgers et al.l 
1196011 and located in the Carina-Sagittarius spiral arm 
(a, 6) = (10'^2 3""58M,-57°45'49") (J20001, {l,b) = 
(284.3°, —0.34°). iMoffat et al.l (jl991ari suggest that Wd2 
contains more than 80 0-type stars. 

Very few young (< 5 Myr) and massive (> 10"'^ Mq) 
star clusters are known in the MW, but they are 
frequently observed in the Magellanic Clouds (e.g., 
iCascoigne fc Kronlll952t iClark et all 1200,51 1 as well as in 
the more distant U niverse^ espec ially in disk and star- 
burst galaxies (e.g.. lHodgelll96l[ ). Therefore, Wd2 is a 
perfect target to study the star formation process and 
feedback of the gas in the presence of massive stars as 
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well as the possible triggering of star formation in the sur- 
rou nding molecular cloud (suggested to occur in R CW 49 
by I Whitney et al.l 12004 iChurchwell et al.l 12004 using 
Spitzer mid-IR images). In fact, man y globular clusters 
are le ss massive than Wd2 (see, e.g., iMisgeld fc Hilkerl 

Honl). 


The Wd2 cluster has been widely discussed in the lit¬ 
erature in the last decade, yet its physical properties 
are poorly known. There is a considerable disagreement 
about the dis t ance o f Wd2, ranging from values of 2.8 kpc 
(jRauw et al.l 120071 lAscenso et al. I2007t iCarraro et al.l 
l2013fl. 4.16 knc (iVargas Alva rez et al .l l2013fl. 5.7 k pc 

(|Piatti et al.lll998ll. 6.4 kncjCarraro fc Munarill2004l to 

even8kDC (iRauw et al.ll2007Ll2011ll . A consensus has not 
been reached on the age of Wd2 either, although there 
seems to be general agreement that the cluster as a whole 
is younger than 3 M yr and that the core may be even 
youn ger than 2 Myr (jAscenso et al.ll20(T7l: ICarraro et al.l 
l2013f l. This young age combined with the cluster half- 
light radius of about 2.4 pc and its total mass make it 
likely that Wd2 i s younger than half of its typical crossing 
time (^5.5 Mvr lGieles fc Portegies Zwartll201ll) . There¬ 
fore, the physical conditions of Wd2 are very close to its 
initial conditions. 

Over the last years a series of spectroscopic ob- 
servations of the m ost massive stars was performed. 
IMoffat et al.l (Il991bll used UBV photometry and low 
resolution spectroscopy in order to classify six 0-stars 
within Wd2 as 06-7 V objects and to obta in a distance 
of d = 7.9tlo while iPiatti et al.l ( 19981 1 reanalyzed a 
subset of O stars from (jMoffat et al.lll99fbll to determine 
a distance of d = 5.7 ± 0.3 kpc. Two snapshot spectra of 
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WR20a were taken and analyzed bv iRauw et al.l (|2004D . 
They found that WR20a is a double-lined binary with a 
preferred orbital period of 3.675±0.030 days with compo¬ 
nent minimum masses of 68.8±3.8 Mq and 70.7±4.0 Mq. 
They su ggest a spectral type of WN6ha for both com¬ 
ponents. [BonanosIer^I] (l2004fl discovered it was also an 
eclipsing binary and refined those values to even more ex¬ 
treme masses of 83.0 ± 5.0 Mq and 82.0 ± 5.0 Mq since 
they also measured the inclination. Spectroscopic ob¬ 
servations of iVargas Alvarez et al.l (j2013fl ; iCarraro et al.l 
(|2013f) of the most massive stars show that Wd2 has a to¬ 
tal to selective extinction in the visual {Rv = Ay/E{B — 
R)) of ab out 3.77-3.85, therefore highe r than the Galac¬ 
tic mean (jMathis fc WallenhorstlllQ^ of Ry = 3.1. 

The proximity to our solar system makes Wd2 a good 
target to resolve and characterize the individual com¬ 
ponents of its stellar content. Still, Wd2 has not been 
extensively studied. We attempted to fill the observa¬ 
tional gap by designing a high resolution multi-band sur¬ 
vey of the Westerlund 2 region with the Hubble Space 
Telescope. In this paper, we present new high reso¬ 
lution photometry, obtained wit h the Advan ced Cam¬ 
era for Surveys (hereafter ACS; lUbedal [201211 and the 
infrared channe l of the Wid e Field Camera 3 (here¬ 
after WFC3 /IR: lDressell2012[) of the Hubble Space Tele¬ 
scope (HST). Combining four wide band filters (F555VF, 
F814IF, F125kF, F160W) we aim to constrain the dis¬ 
tance and age of the cluster and to characterize its low- 
mass stellar content. Combining two narrow band filters 
{F658N, F128N) we can d erive a pixel-to-pixel redden¬ 
ing map (jPang et al.ll201l!) with the Ha to Pa,5 ratio to 
get the F{B — V) color excess of the gas content in Wd2. 
This map helps us to deredden the stellar content while 
taking into account differential reddening. 

In Section[2]we describe the HST observations of Wd2, 
and the data reduction, bas ed on the HST pipelines 
(|DolDhinll2000l : [lYuchtei]l201^ . In Section |3] we describe 
the creation of our photometric catalog as well as the 
cleaning process from artifacts and spurious objects. In 
Section [1] we create the high resolution two-dimensional 
map of the color excess. In Section [5] we describe the 
transformation of the gas excess to a stellar excess. In 
Section [6] we describe the morphology and limitations of 
our color-magnitude diagram. In Section [7] we determine 
the physical parameters of Westerlund 2 such as distance, 
age, and the total-to-selective extinction towards Wd2. 
In Section [H we present the unusual morphology of Wd2 
that suggests that Wd2 might actually be a double clus¬ 
ter. In Section[3we summarize findings and conclusions. 

2. HST OBSERVATIONS AND DATA REDUCTION 
2.1. Observations 

Our observations of Wd2 were p erformed with the HST 
during Cycle 20 using the ACS (lUbedal 1201211 an d the 
Infrared Channel of the WFC3/IR ([P resse f lmi . Six 
orbits were granted and the science images were taken 
on 2013 September 2-8 (ID: 13038, PI: A. Nota). 

We acquired data in six different filters: four wide band 
filters and two narrow band filters. We imaged Wd2 
with ACS (resolution 0.049"/px) through the F555IU 
and F81AW filters, corresponding to the V and /-bands, 
as well as the F658N hlter, centered on the Ha line. 
We also observed the cluster with WFC3/IR (resolution 


0.135"/px) and the F125IU and F160IU filters matching 
the J and //-band, as well as the F128N filter corre¬ 
sponding to the Pa/3 line. In addition, short exposures 
were acquired in the /"555IF and F81AW bands (here¬ 
after denoted as F555IU(s) and F814IU(s)) to reduce 
saturation due to 0-type stars in the very bright cluster 
center and across the field of view. 

We used one ACS pointing to cover Wd2 and a dither 
pattern with four positions for the long-exposures (350 s 
each) in the broad band filters, subsequently replaced 
by a two-position dithering for the F658N filter (for a 
total exposure time of 1400 s). The short exposures were 
only observed at one position, each with an exposure 
time of 3 s. For the WFC3/IR observations mosaics of 
2x2 tiles were required since the WFC3/IR has a smaller 
held-of-view (FOV) than the ACS (FOVacs = 202 x 
202 arcsec and FOVwfcs = 123 x 136 arcsec). Each 
wide-band tile consists of four dither positions with a 
total exposure time of 947 s (3x250 s and 1x200 s). The 
narrow band tiles consist of three dither positions with 
a total exposure time of 748 s (250 s each). The log 
of the observations is reported in Tab. [T] together with 
the filters and exposure times used. An overview of the 
spatial positions of all images as well as their alignment 
is shown in Fig. [1] 



RA 


Fig. 1.— The field of view for the different exposures taken 
with the WFC3/IR and ACS filters. The blue and black fields 
together represent the F125W and FldOW WFC3/IR exposures, 
while the blue fields alone represent the F128N WFC3/IR dithered 
images. The yellow, green, and red fields together represent the 
ACS ^8141^ and Fb5bW long exposures, while the yellow fields 
represent the F658N dither positions. The red field alone repre¬ 
sents the short exposures for ACS F814W and FbbbW. 


2.2. Data reduction 

The obser ved data were first process ed with the inter¬ 
nal pipeline (jIJbedall2012l : lDresse]|[2012h (from now on re¬ 
ferred to as raw data), which delivers pipeline flat-fielded 
files (fit.fits for WFC3) or charge-transfer efficiency 
(CTE) corrected images (flc.fits for ACS). 

The reduction of the raw data consists of mu ltiple 
steps performed by AstroDrizzle (|Fruchtei] 120101) . the 
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TABLE 1 

Overview of the observations 


Dataset 

Filter 

Camera 

RA(J2000) 

Dec(J2000) Exp. Time [s] # Exp. 

Date 

IBYAOlOlO 

F555W 

ACS 

10:24:08.989 

-57:46:20.95 

1400 

4 

2013-09-02 

IBYAOlOlO 

F555IF(s) 

ACS 

10:24:08.989 

-57:46:20.95 

3 

1 

2013-09-02 

JBYA01020 

Fe58N 

ACS 

10:24:03.756 

-57:45:33.02 

1400 

2 

2013-09-02 

JBYA01030 

F81AW 

ACS 

10:24:08.989 

-57:46:20.95 

1400 

4 

2013-09-02 

JBYA01030 

F814VF(s) 

ACS 

10:24:08.989 

-57:46:20.95 

3 

4 

2013-09-02 

IBYA08010 

Fl2bW 

WFC3/IR 

10:24:07.294 

-57:44:05.26 

947 

4 

2013-09-02 

IBYA09010 

F125W 

WFC3/IR 

10:24:16.116 

-57:45:48.56 

947 

4 

2013-09-02 

IBYAlOOlO 

F125W 

WFC3/IR 

10:24:55.869 

-57:45:07.69 

947 

4 

2013-09-03 

IBYAllOlO 

F125W 

WFC3/IR 

10:24:04.686 

-57:46:51.05 

947 

4 

2013-09-08 

IBYA08030 

F128N 

WFC3/IR 

10:24:07.294 

-57:44:05.26 

748 

3 

2013-09-02 

IBYA09030 

F128N 

WFC3/IR 

10:24:16.116 

-57:45:48.56 

748 

3 

2013-09-02 

IBYA10030 

F128N 

WFC3/IR 

10:24:55.869 

-57:45:07.69 

748 

3 

2013-09-03 

IBYA11030 

F128N 

WFC3/IR 

10:24:04.686 

-57:46:51.05 

748 

3 

2013-09-08 

IBYA08020 

Fl&OW 

WFC3/IR 

10:24:07.294 

-57:44:05.26 

947 

4 

2013-09-02 

IBYA09020 

Fimw 

WFC3/IR 

10:24:16.116 

-57:45:48.56 

947 

4 

2013-09-02 

IBYA10020 

Fl&OW 

WFC3/IR 

10:24:55.869 

-57:45:07.69 

947 

4 

2013-09-03 

IBYA11020 

Fimw 

WFC3/IR 

10:24:04.686 

-57:46:51.05 

947 

4 

2013-09-08 


Note. — This table summarizes the observations of Wd2 taken with the ACS and WFC3/IR cameras 
aboard HST. Given are the identifiers of the datasets, the filter, camera, as well as the central coordinates of 
each dither position. The exposure time is the total exposure time of all dither positions combined. We also 
give the observation date for each image. 


follow-up tool of MultiDrizzle in the HST pipeline. The 
major goal of AstroDrizzle is to reconstruct and re¬ 
cover the spatial information of a set of images while 
preserving the signal-to-noise ratio (SNR). In Astro- 
Drizzle this is done by the method formally known as 
variable-pixel linear reconstruction (commonly known 
as Dri zzle) developed by An dy Fruchter and Richard 
Hook (|Fruchter fc HookI 1199711 . A detailed description 
of the package can be found in the manual of DrizzlePac 
(http://drizzlepac.stsci.edu). 

AstroDrizzle performs its tasks on pipeline flat-fielded 
files or CTE corrected images. The general concept of 
AstroDrizzle is to partly recover the high frequency infor¬ 
mation smeared out by the detector pixel response func¬ 
tion by combining the sub-pixel dithered images. This 
algorithm gives us the chance of recovering spatial infor¬ 
mation of the under-sampled WFC3/IR images to recon¬ 
struct a final combined image with a higher resolution. 
In the case of multiple dither-positions the final com¬ 
bined (drizzled) image is free of cosmic rays and detector 
defects. At the end all images are aligned to the same ref¬ 
erence frame, provided by the world coordinate system 
(WCS) information in the header of one of the images 
(we chose the F814W filter as reference). 

Since the instruments (ACS and WFC3) and the 
pipeline data products are substantially different we treat 
the related images separately. 

2.2.1. The ACS data 

During the whole process the short exposures for the 
F'o'b'oW and the FSldkF filters are treated individually 
and are labeled as E555IF(s) and F814VF(s), since As¬ 
troDrizzle cannot effectively combine images of very dif¬ 
ferent exposure times. 

For all ACS images we used a pxfrac=0.9 for the driz¬ 
zling process meaning that each pixel is shrunk by 10% 
in size before being mapped onto the output grid. This 
value was chosen so that the size of the point-spread- 
function (PSF) of the drizzled image reaches a minimum 


while mapping all pixels properly to the output grid (see 
Sect. 12.31) . The output pixel grid has the same pixel size 
as the original image (0.049 arcsec/pixel). Since there 
is only one dither position available for the short expo¬ 
sures (see Fig. [T]), the cosmic-ray and detector defect re¬ 
moval could not be accomplished. Nevertheless we driz¬ 
zled those images to correct them for geometric distortion 
and to align them properly to the reference frame. The 
short exposures are only taken to recover the photometry 
of sources that are saturated in the long exposures. 

2.2.2. The WFC3/IR data 

The PSF of the IR data is under-sampled. Therefore, 
we reduced the pixel size of the output grid to recover 
as much spatial information as possible for the deter¬ 
mination of the optimal pixel size. If the pixel size in 
combination with the pixel fraction is shrunk too much, 
there is not enough information left for mapping the 
pixels properly onto the output grid without generat¬ 
ing ’’holes” of empty pixels. Taking into account the 
weighted maps while representing the contribution of 
each pixel to the final image and the final full width at 
half maximum (FWHM) of the PSFs of some represen¬ 
tative stars, we chose a final pixel size of 0.098 arcsec/px 
with a pxfrac=0.8. The data would have allowed us to 
reduce the pixel size a little more, but we found it con¬ 
venient to have a pixel size twice that of the ACS. In 
addition, we mosaicked the four pointings into a single 
image. 

In the end we had five mosaics for the ACS data set, 
one for the narrow-band and two (short and long expo¬ 
sures) for the wide-band filters. For the WFC3/IR filters 
we had three mosaics in total, one for each filter. 

2.3. Quality check of the final images 

In order to check the quality of the drizzled ACS im¬ 
ages, we tested the stability of the weighted maps of the 
final products. In order to determine the variations in 
the weighted maps, 28 boxes (20 for the F658N filter) 
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TABLE 2 

Summary of the weighted maps of the AstroDrizzle 

PIPELINE 


Filter 

# Boxes 

Boxsize [px] 

median 

(T 

(T 

median 

F5551F 

28 

100 X 100 

1318 

147.7 

0.1147 

F555W 

16 

250 X 250 

1325 

142.2 

0.1099 

Fe58N 

20 

100 X 100 

1336 

208.5 

0.1618 

F658N 

16 

250 X 250 

1339 

206.7 

0.1599 

F814:W 

28 

100 X 100 

1317 

147.3 

0.1145 

F814W 

16 

250 X 250 

1322 

147.7 

0.1137 


Note. — In this table we give a summary of the tests for the 
weighted maps. The first column shows the corresponding filter. 
Column 2 gives the number of boxes per image. Column 3 lists 
the boxsizes for calculating the median and the standard devia¬ 
tion of the pixel values. Column 4 and 5 show the median and 
the standard deviation of the pixel values, respectively. Column 
6 lists the ratio of the standard deviation and the median. 

of 100 X 100 pixel and 16 boxes of 250 x 250 pixel were 
selected across the image. In each box the standard de¬ 
viation of the median was determined using the IRAFQ 
task IMEXAMINE. This ratio should be below 0.2 following 
the DrizzlePac manual, meaning that the deviation of the 
weights of the different pixels should be less than 20%. 
This test could not be performed for the IR data due to 
the different readout method. The results are shown in 
Table m 

In addition, we determined the FWHM of a sample of 
25 carefully selected stars across each field with the IRAF 
task IMEXAMINE. We chose medium bright, isolated, non- 
saturated stars so that the PSFs are not contaminated by 
close neighbors or remaining detector artifacts. For the 
PSF determination, all ACS images where divided into 
four quadrants. The mean FWHM of the 25 selected 
stars in each filter is listed in Tab. [31 The FWH M range 
over the whole chip is given in the ACS manual (lUbedal 
IWll with the FWHM of F555M = 2.12-2.49 pixel as 
reference. The FWHM is quite constant between the dif¬ 
ferent ACS filters and even the WFC3/IR filters. Our de¬ 
termined values are within the PSF range in the F555M 
filter, meaning that quality and, therefore, the drizzle 
parameters are fine. 

In Fig. [2] we show a color composite image of the Wd2 
region, where the F125IF, F814IF, and F555IF filters 
are represented in red, green, and blue, respectively. 

3. THE CREATION OF THE PHOTOMETRIC 
CATALOG 

The source detection and photometry was performed 
with the DOLPHO40 package including the extensions 
for ACS and WFC3/IR written by Andrew Dolphin. 
DOL PHOT is a modified version of HSTphot (iDolohinl 

3.1. Source detection and photometry 

For the source detection and photometry we adopted 
individual parameters for each filter. We considered the 

^ IRAF is contributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

® http://americano.dolphinsim.com/dolphot/ 


size of the PSFs, the different saturation limits, and the 
nebulosity conditions (due to the very different level of 
gas and stellar density). The main effort was to detect 
all sources and yet to avoid as much contamination by 
spurious objects as possible in order to make the clean¬ 
ing procedure later-on easier. Spurious detections may 
happen in particular in regions with varying background 
brightness due to the varying gas and dust density in the 
H II region and stellar crowding. 

DOLPHOT uses the drizzled images for the source 
detection and works on images (Acs and fits) with 
geometric-distortion-corrected headers (performed by 
AstroDrizzle) (from now on called processed images). To 
perform the source detection several steps have to be ap¬ 
plied in a given order: acsmask (wfcSmask) to mask the 
processed images for detector artifacts and other defects 
using the weighted images, splitgroups to create single 
chip and single layered fits files, and calcsky to create a 
skymap for all processed images. 

After those preparatory steps the actual dolphot rou¬ 
tine was used to create a point-source catalog using the 
PSF and pixel area maps from an analytic model of ACS 
and WFC3/IR which is provided on the DOLPHOT web¬ 
page. The photometric run was performed in one step 
using the reference frame (F814IF) to obtain the com¬ 
mon astrometry and in order to improve the photometric 
accuracy. The output of DOLPHOT is one catalog for 
each filter and exposure, so in total 8 catalogs. 

3.2. Improving the quality of the Photometric Catalogs 

Initially, the photometric catalogs contain a number of 
spurious objects, which must be removed. We ’’cleaned” 
the photometric catalogs by applying appropriately cho¬ 
sen trimming parameters for each of the filters. The 
trimming values for the WFC3/IR and ACS instruments 
are given in the DOLPHOT manual. In order to ensure 
that as many spurious objects as possible were deleted, 
still preserving all ’real’ objects, especially in the nebu¬ 
lar regions and the dense cluster center we had to find 
the optimal parameter set for DOLPHOT. These values 
were adjusted by visual inspection of the observations of 
the removed and the kept objects. 

3.2.1. Object Type 

The objects are classified as different types (OT) by 
DOLPHOT: good star (1), star too faint for PSF deter¬ 
mination (2), elongated object (3), object too ’’sharp” 
(4) (see Sect. 13.2.311 . extended object (5). We only kept 
the objects with OT=l. For the WFC3/IR filters we did 
not select by object type in order not to lose any objects. 
The distribution of the OTs for the ACS filters can be 
seen in Fig|3l 

3.2.2. Magnitude error 

We restricted the catalogs to contain only those stars 
with magnitude errors above a given threshold, which 
depends on the filter. The cuts for the different filters 
are listed in column 2 of Tab. HI 

3.2.3. Sharpness 

The ’’sharpness” parameter allows us to assess how 
well the PSF of a star is fitted. For a perfectly-fitted 
star the value is zero, positive for a point-source that is 
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TABLE 3 

The mean PSF FWHM in the different filters 


Filter 

FWHM [px] 

<^FWHM [px] 

FWHM [arcsec] 

cpwHM [arcsec] 

F555W 

2.154 

0.1044 

0.1055 

0.0125 

F658N 

2.034 

0.3848 

0.0997 

0.0187 

F81AW 

2.011 

0.2304 

0.0986 

0.0113 

^12514" 

2.188 

0.2145 

0.2144 

0.0210 

^128^ 

2.410 

0.2362 

0.2362 

0.0231 

FIQOW 

2.242 

0.2197 

0.2197 

0.0215 


Note. — The mean FWHM of the PSF of 25 selected stars in each filter. 
The first column shows the corresponding filter. Column 2 and 3 show the 
FWHM and the standard deviation in pixels and column 4 and 5 show the 
FWHM and the standard deviation in arcsec. 


TABLE 4 

The thresholds for cleaning the catalogs 


Filter 

u [mag] 

Sharpness 

Crowding 

F555VF 

0.5 

±0.2 

0.10 

F555W(s) 

0.5 

±0.2 

2.50 

F658N 

0.5 

±0.4 

0.10 

F81AW 

0.5 

±0.4 

3.50 

F81AW{s) 

0.5 

±0.2 

2.50 

F125W 

0.2 

±0.5 

0.15 

F128N 

0.2 

±0.5 

0.20 

Fimw 

0.1 


0.20 


Note. — In this table we give the different 
thresholds for cleaning the photometric catalogs 
from spurious objects. The columns 2-4 refer to 
the quality parameters for point-source photome¬ 
try as provided DOLPHOT. Objects with values 
equal to or smaller than the listed thresholds were 
retained. 

too peaked as compared to a typical stellar PSF such 
as a cosmic ray, and negative for a point source that is 
too broad in its light profile and that might be a back¬ 
ground galaxy. Typical values range from —0.3 to -1-0.3 
for un uncrowded region. The values chosen as our se¬ 
lection criteria are given in column 3 of Tab. |4l For the 
F128N filter a cut was impossible without losing ’’real” 
objects. The distribution of the sharpness as a function 
of magnitude can be seen in Fig. 01 

3.2.4. Crowding 

The crowding is a measure of the flux contamination 
by surrounding objects. For an isolated star the value 
is zero. Crowding effects become especially important in 
very densely populated regions or in areas around very 
bright saturated stars since the crowding parameter gives 
a value in magnitudes by how much brighter a star would 
have been measured if nearby stars would not have been 
fit simultaneously. Our values chosen to select only those 
stars with a well-defined PSF and very low contamination 
can be seen in column 4 of Tab. 01 This selection is 
visualized in Fig. IS] 

3.3. Merging the different filters 

For the creation of the final catalog we merged the dif¬ 
ferent filter catalogs. We used the pixel coordinates of the 
resampled images and we established a maximum radial 
distance of 1.5 pixels (0.0735 arcsec/pixel) for merging 
the different filters. This means that the position of an 


object in two different filters is not allowed to radially 
deviate by more than 1.5 pixels from each other. 

First, we combined the long and the short exposures 
for the F5551F and the F8141F band. Then, we used the 
F814W catalog as basis, and subsequently added one fil¬ 
ter after another. At the end, we adopted the World Co¬ 
ordinate System (WCS), in J2000 coordinates, which is 
provided by the processed long exposures in the FSIAW 
filter as the common coordinate system for all exposures. 

In order to merge the long and short exposure cata¬ 
logs this procedure was repeated with a maximum radial 
distance of 1 pixel, between the same object in the short 
and long exposures. Objects from the short exposures 
were only taken into account for magnitudes brighter 
than 20 mag. For fainter objects, the long exposures 
were not saturated and offered a higher signal-to-noise 
(SNR) than the short exposures. 

For all objects available in both exposures 
(''t-F 555VF=246 and npsi4w=(i73, where n is the 
number of objects) we chose to adopt the point-source 
photometry with the smaller photometric error which 
always came, as expected, from the long exposures. 

After all filters had been merged into one catalog, all 
objects that were not detected in at least two filters were 
deleted. It is unlikely that a star is only detected in one 
filter. This is not true for the F160W band, however, 
since it can happen that a ve ry young object emits only 
far to near-infrared radiation (iLada fc Wilkindfl98411 and 
therefore is only detected in the i7-band. In order not 
to lose those (potentially interesting) objects we decided 
not to remove the F160IF-only detections. 

At the end we have a catalog consisting of objects de¬ 
tected at least in two filters (except for the F160IF filter) 
with a common coordinate system in pixel coordinates, 
as well as a WCS based on the F8I4kF images. 

3.4. Properties of the photometric catalog 

The final catalog consists of 17121 objects in total. 
For an overview of the magnitude range as well as of the 
depth of each band, we refer to Fig01 The arrows mark 
the limiting magnitudes up to where 90% of all sources 
are located. A summary of the number of point-source 
detections in the combined catalog is given in Tab. 0] 

The photometric errors are, of course, below the cuts 
set in the cleaning process. As one can see in Fig. 0 
most sources have much smaller errors than these limits. 
The limiting errors and resulting magnitudes comprise 
90% of all detected point sources in a given filter when 
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Fig. 2.— Color composite image of the HST ACS and WFC3/IR data of Wd2, including the ^12514^ (red), ^81414^ (green), and F5554F 
(blue) filters. North is up. East to the left. The FOV is ~4 arcminx4 arcmin. This image was chosen to be the official Hubble 25th 
anniversary imag43. Credit: NASA, ESA, the Hubble Heritage Team (STScI/AURA), A. Nota (ESA/STScI), and the Westerlund 2 
Science Team. 

®http://hubblesite.org/newscenter/archive/releases/2015/12/image/a/ 


no cuts are imposed, and are given in Tab. [S] 

3.5. The photometric improvement 

The photometric data of our survey are the deepest 
data obtained for Wd2 so far. The Wid e-Field Plane¬ 
tary Camera 2 (WFPC2) (|Gonzagall201C)D HST data of 
IVargas Alvarez et al.l (I2013D reach a photometric uncer¬ 
tainty level of 0.2 mag at a depth of ^ 21.5 mag and 


~ 21.0 mag in the F5551T and F8141T filters, respec¬ 
tively. At the same level of photometric uncertainty 
we go down to ~ 27 mag in the F5551T filter and to 
~ 26 mag in the F8141T filter, implying that our data 
are > 5 ma g deeper in th e optic al. The J and H pho¬ 
tometry bv lAscenso et al.l (1200711 has a limiting magni¬ 
tude of ^ 20 mag and ~ 18.5 mag at a photometric un¬ 
certainty limit of 0.1 mag, respectively, where our HST 
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Fig. 3. — The total number of stars for different object types 
per ACS filter in the catalogs created by DOLPHOT. We decided 
to keep only those objects with OT=l to reduce the number of 
spurious detections. 


TABLE 5 

Properties of the final catalog 


Filter 

Total Number 

mag-LiMiT [mag] 

cr-limit [mag] 

F555VF 

3 922 

26.7 

0.160 

F658Af 

3 556 

23.8 

0.185 

F814VF 

9 281 

25.3 

0.093 

F125IF 

10 981 

23.1 

0.035 

F128Af 

8 955 

20.8 

0.095 

F160W 

15 355 

22.6 

0.054 


Note. — In this table we give the main properties of the final 
photometric catalog. In column 2 we show the total number 
of all objects per filter. Column 3 and 4 give the limits in 
magnitudes and the photometric error below which 90% of all 
sources are located. In total, 17121 objects were detected. 2 236 
point sources had a detection found in all filters. 


photometry in the F\2'oW and F160W filter goes down 
to ^ 23.5 mag and ~ 23 mag at the same uncertainty 
levels. Also in the near infrared (NIR) filters our data 
go ~ 3-5 mag deeper than the photometric catalog of 
IVargas Alvarez et al.l (|2013ll . Here it is worth reminding 
the reader that even when including the 3 s short expo¬ 
sures of the ACS broad-band filters we miss 16 bright 
objects (see Subsect. OD . Those objects are too bright 
or blended to be accurately measured by DOLPHOT. 

The spatial resolution of our dataset is also unprece¬ 
dented. Especially due to the sub pixel drizzling we are 
able to reach a resolution of 0.049 arcsec pixeH^ in the 
optical and 0.098 arcsec pixel“^ in the IR. This is al- 
m ost 50% highe r than the pixel resolution of the data 
of lAscenso et al.l (|2007ll (0.144-0.288 arcsec pixeD^) and 
m ore than eight t imes higher than the optical resolution 
of iCarraro et al.l (|2013ll (0.435 arcsec pixel'^). 

These deep high resolution data in 6 filters give us the 
unique opportunity to study the properties of Wd2 at an 
unprecedented level of accuracy and depth. They reveal 
the young low-mass population (see Sect. El), give us 
new information on the spatial distribution of the cluster 
members (see Sect. 0), make it possible to create for the 
first time a high resolution pixel-to-pixel map of the color 
excess in order to correct for differential reddening (see 


Sect. [3]) and to improve the cluster age determination. 

4. A HIGH-RESOLUTION TWO-DIMENSIONAL 
MAP OP THE COLOR-EXCESS 

Wd2 is located in the Carina-Sagittarius arm of the 
MW, and thus the local extinction due to interstellar 
dust is expected to be high. A high resolution, spatially 
resolved characterization of the extinction in the region 
is key to a precise determination of age and distance 
for the Wd2 cluster, and to investigate the existence of 
an age spread. F ollowing the methodology outlined by 
IPang et all (IWl for NGC3603, we used a combination 
of Ha and Pa/3 images {F658N and F128N) to derive a 
high-resolution, pixel-to-pixel map of the color excess of 
the gas (g) E{B — V)g of the Wd2 region. 

4.1. Ha and Pajd emission 

In order to construct the pixel-to-pixel color excess 
map, we exploited the fact that radiation emitted at 
shorter wavelengths is absorbed more easily by dust than 
radiation emitted at longer wavelengths. Therefore, the 
interstellar dust attenuates the light from a given source 
emitted in the Ha line (rest-frame wavelength 6 563 A) 
more than the light emitted in the Pa^ line (rest-frame 
wavelength 12 802 A). Consequently, the observed flux 
ratio of Ha/Pa/3 is always smaller than the theoretical 
(reddening free) ratio computed for the same conditions 
o f electron densi t y and electron temperature. 

iCalzetti et al.l l| 199611 provides a relation between the 
color-excess E{B — V)g of the interstellar gas and the 
observed {Robs) and theoretical (Rmt) flux ratio Ha/Pa/3 
with an assumed extinction law, represented by the total 
to selective extinction k(A) = A{X)/E{B — V) for a given 
wavelength A: 


E{B - V) 


log {Robs IRint ) 

0.4 [At(Ha) — K(Pa/3)]' 


( 1 ) 


In order to use Equation El we need to get the pixel- 
to-pixel Ha/Pa/3 flux ratio from our observations. We 
applied the following steps in order to properly remove 
the stellar emission: 


• First, we repeated the drizzling process (see Sec. 
ED for the ACS images adopting the same pixel size 
as for the WFC3/IR images (0.098 arcsec pixel“^). 
A visual inspection ensured that the images where 
properly aligned. 

• Due to the use of different instruments covering 

the whole wavelength range we need all of our im¬ 
ages to have the same PSFs. Therefore, we used 40 
stars that were non-saturated and without a very 
close neighbor in each of the six filters in order to 
determine the PSF. The PSF was largest in the 
E128N image with a FWHM of 2.26 ± 0.312 pix¬ 
els. Hence we degraded all the other images to the 
same resolution by convolving them with a Gaus¬ 
sian whose dispersion (cr^) is the difference between 
the FI28A^ filter (o'|.;^ 28 Af) filter 

in question. 


• We calibrated all images in flux by multiplying 
them by their specific PHOTFLAM value given in 
the respective header. 
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Fig. 4.— The distribution of sharpness, which is a measure of how peaked the source is compared to a stellar PSF as function of magnitude 
for each filter. The black dots represent the data retained. 


• The continuum flux centered at the Ha and Pa/3 
wavelength was derived by linear extrapolation of 
the fluxes in the F555tP and F8141P {Fl2bW 
and F160W) image pixel-by-pixel at the F658N 
{F128N) wavelength. With such an extrapolation 
we take into account the slope of the continuum 
emission of the stars in order to better remove them 
from the final, pure line-emission images. 

• The interpolated continuum image was then sub¬ 
tracted for properly removing the stars from the 
F658A^ (F128iV) image. Afterwards each pixel 
value was multiplied by the respective bandwidth 
resulting in the pure line emission image as well as 
the Robs flux ratio. 

• For the theoretical flux ratio of Ha/Pa/3, Rmt, 
we derived a value of 17.546 with an assumed 
electron temperature of 10 000 K and an elec¬ 
tron density of 100 cm“^ (jOsterbrock fc Bochkarev! 
119891 1. We adopted a normal total-to-selective ex- 
ti nction of Ry = 3.1. From the extinction law 
of iCardelli et al.l (jl989H we derived KfHq) = 2 437 
and K(Pa/3) = 0.8035. With ll^zpatrickl (1199911 we 
derived ^(Ha) = 2.346 and K(Pa/3) = 0.7621. 

• Finally we were able to apply Equation [T] to con¬ 
struct the pixel-to-pixel map of E{B — V)g shown 


in Figure [H 

4.2. The features of the color excess map 

The two-dimensional color-excess map (see Fig. [8]) is 
created with a spatial resolution (0.098 arcsec pixel~^) 
that has never been reached before for this region. The 
very inhomogeneous distribution of the dust extinction 
shows that differential reddening must be taken into ac¬ 
count before any conclusion can be reached on the phys¬ 
ical parameters of Wd2. We will use this map to get 
a dereddened cluster population, where for each of the 
stars the individual reddening is estimated locally (see 
SectIZl). 

In Fig. [9l the reddening map is shown on the left in 
black and white over-plotted with the selected cluster 
members (see Sect. [Zl) in green, and the contours rep¬ 
resenting E{B — V)g. On the right we show the RGB 
image of the same spatial extent over-plotted with the 
same E{B — V)g contours. This image emphasizes the 
location of the stellar population of Wd2 with respect to 
the highly variable color excess on small scales. Here it 
becomes clear that this map will improve the correction 
for differential reddening of each star individually. 

The color excess map is morphologically very inho¬ 
mogeneous. In the center where the cluster is located 
(marked by a cross in Fig. |8]), E{B — V)g has a value of 
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Fig. 5.— The crowding as a function of magnitude for each filter. The black dots represent the data retained. 


^1.8 mag. Here we also see most of the very bright stars 
in the color composite image (see Fig. El). Also in the 
very central region around the cluster there still seems 
to be a concentration of gas and dust, although we can¬ 
not say whether these components are mainly located 
in front of the cluster as seen along our line of sight, 
or in the cluster. A ’’bridge” with a slightly increased 
E[B — V)g w 1.9 mag extends from the cluster center 
to the North-East. In Fig. El one can see that there are 
no massive luminous stars visible within this ’’bridge”. 
This ’’bridge” is not visible in the RGB image. Towards 
the South and the East where a higher concentration of 
dust and gas is observed, the color excess increases to 
values larger than 2.0 mag. In very dense regions it even 
can exceed 2.5 mag. The color excess map follows very 
well the morphology of the gas clouds showing a direct 
connection between the color excess E{B — V)g and the 
cloud structure. As expected the color excess is higher 
where the clouds are denser. 

The median reddening of the pixels across the redden¬ 
ing map is E{B — V)g = 1.87 mag. The vast majority of 
all pixels have E{B — V)g between 1-3 mag. Only 10 742 
pixels (0.23%) have E{B — V)g > 3.0 mag and only 1 018 
pixels (0.02%) exceed 5.0 mag. The distribution of the 
nebular color excess values can be seen in Fig. (TU] 

As shown in Fig. |8] the range of E{B — V)g is quite 
large. We carried out an analysis of the radial distribu¬ 


tion of the color excess. We divided the color excess map 
into four quadrants (NW, NE, SW, SE) (see Fig. [8]), with 
the origin at the center of the cluster (cross in Fig. HI). 
In Fig. ITlIone can see the azimuthal dependence of the 
color excess. The dashed-dotted line shows the average 
of all quadrants. To the South the color excess increases 
almost constantly from ~ 1.86 mag to 2.15 mag in the 
South-West and ^2.0 mag in the South-East. To the 
North, E{B — V)g decreases to ^ 1.7 mag. 

4.3. Uncertainties of the color excess map 

Before applying the reddening map to our photomet¬ 
ric measurements, we needed to assess the accuracy of 
the map itself. We did that by assuming that the elec¬ 
tron counts on the detector follow a Poisson distribution. 
We started with a Poisson error of the counts per pixel 
and performed a proper error propagation throughout 
our whole computation, beginning with the continuum 
interpolation up to the final E(B — V)g from Equation 

[H 

In Eig. [12] we show the distribution of (Je{b-v) with 
a median of 0.032 mag. 90% of all pixels have an er¬ 
ror smaller than 0.057 mag in the inferred nebular color 
excess. 

As mentioned before, the reddening map depends on 
the adopted extinction law. For all further analyses we 
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Fig. 6. — The magnitude distribution of the detections contained in the final Wd2 catalog. The plateau at the faint-end detections in 
the F1601T filter is due to the inclusion of all single-band detections. The arrows mark the limiting magnitudes up to where 90% of all 
sources are located (see Tab. 0. 
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Fig. 7. — The photometric errors in each filter. The horizontal 
lines show the limits below which 90% of all sources are located in 
the deep exposures (see Tab. 0. 

used the extinction law of iCardelli et alJ (|1989L hereafter 
denoted with the subscript C). For the sake of compari¬ 


son, we als o created the color excess map with the extinc¬ 
tion law of lFitzpatricld (|1999l hereafter denoted with the 
subscript FP). Both extinction laws are almost identical 
in the optical regime and so are the results. The median 
of the difference between E{B — V)fp and E{B — V)c 
is 0.05 ± 0.018 mag. 

A non-negligible effect is the contamination of the 
wide-band filters with emission lines (e.g.: O III, S III, 
H/3 ). Due to a lack of spectra of the gas and the large 
differences in published line emissions between different 
H II regions, we are unable to estimate the amount of 
contamination due to those lines. Compared to the Pa/3 
line emission we expect this contamination in the wide¬ 
band filters to be small. 

In order to visualize the possible contamination of the 
color excess map by nebular line emission we show in 
Fig. [T3]the throughputs (top panel) of our filters together 
with the simulated spectral lines (bottom panel) of a typ¬ 
ical H II region usin g the H II Regions Li brar^ based on 
the code CLOUDY (iPanuzzo et al.l[200^ . For the model 
parameters we adopted the values (rtn ^ 100 cm“^, 
Qo = 4.2 • 10®^ s“^) of iPelleg-rini et al.1 ()2nilL and ref¬ 
erences therein) for the massive star cluster R136, since 
its content of massive stars is comparable to Wd2. Since 

® http://pasquale.panuzzo.free.fr/hii/ 
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Fig. 8. — The pixel-to-pixel map of the color excess of the gas E{B — V)g, covering the area where all six filters are available (see 

a [T|l. The cross marks the center of Wd2, as given in Simbad. The green lines indicate the separation in four quadrants (see Sect. 

. E{B — F)g,inin = 0.002 mag and E{B — F)p,max = 11-7 mag. The colorbar shows the color excess in magnitudes. All pixels with 
E(B — V)g > 5 are shown in white. The circular areas near and in the cluster center are caused by the substitution of the saturated pixels 
by the median of the surrounding area. The two slightly inclined short stripes (one on the left, one on the right side) are caused by the 
gap between the two ACS chips (see Fig. 0. They remain at locations in the combined image where only one dither position is available. 
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Fig. 9.— Left: The pixel-to-pixel map of the color excess E{B — V)g covering the area where all six filters are available (see Fig. [T}. 
The green points are the stars in RCW 49 (see Sect. [Til, and the contours represent E{B — V)g in the range of 0 to 2.86 in 9 linear steps. 
Right: RGB composite image of the F125W (red), F8141F (green), and FbbbW (blue) filter. 



Fig. 10.— The distribution of the color excess E(B — V)g with 
a median of 1.87 mag. Only 10 742 pixels (0.23%) have a value of 
E{B — V)g > 3.0 mag and only 1 018 pixels (0.02%) exceed 5.0 mag 

we are only qualitatively interested in the possible gas 
emission lines the adopted values fulfill our needs. 

The Pa/3 line is the only line for which we have informa¬ 
tion in our survey. Ha is located outside the broadband 
filters and, therefore, not contaminating the continuum. 
Assuming the flux density to be the same in F\2bW and 
F128A^, we can derive the contribution of Pa/3 via the 
flux ratio FVlbW / F\2iN . Both filters were multiplied 
by their respective bandwidth and PHOTFLAM value. 
We estimate the contamination of the continuum by the 
Pa/3 to be ~ 22%. 

Taking into account results from iPang et al.l (120111) we 
estimate the contamination of the other possible emission 



radial distance [arcsec] 

Fig. 11.— The averaged radial distribution of the color excess 
E[B — V)g in the four quadrants around the cluster center (cross 
in Fig. EJi including the average over all directions (dashed-dotted 
line). 

lines in the other filters to be about 10-15% which could 
increase the overall color excess by ^0.1-0.15 mag. 

Finally, the E{B — V)g excess is also affected by the 
electron density and temperature used to calculate the 
Rint ratio. We calculated E{B — V)g for different elec¬ 
tron densities and temperatures. For an increase of 
the electron density by a factor of two, E{B — V)g de¬ 
creases by 0.00124 mag. For the extreme cases of Ne = 
10 000(100) cm-3 and = 5 000(20 000) K, E{B - V)g 
changes by -|-0.028(-0.034) mag. 

In summary, we can say that the uncertainties in elec¬ 
tron density and temperature are introducing a small 
spread, which is below the 90% limit of the photomet¬ 
ric uncertainty of the different filters (see Fig. [ 21 ). The 
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Fig. 12.— The distribution of o-E(B—V)g with a binsize of 
0.0075 mag. The median value is 0.032 mag. The number is nor¬ 
malized to the total number of pixels. 
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Fig. 13.— Top: Shown are the throughputs of the six used 
filters, derived with Synphot. The F128A'^ passband is located 
within the wavelength range covered by the F12bW filter. There¬ 
fore, the Pa/3 emission line contributes to the 1^12514^ continuum. 
Bottom: Synthetic spectral lines from the H II Regions Library 
UPanuzzo et al.l2003ll show the emission features contaminating the 
continuum filters. The Ho line is truncated for plotting reasons. 

unknown line contamination, on the other hand, intro¬ 
duces a possible increase of the color excess of the ga,s 
E{B — V)g of ~0.1-0.15 mag. This shift can be cali¬ 
brated out when transforming E{B — V)g to the stellar 
color excess E{B — V)* (see eq. [S|). 

4.4. Limitations of the color excess map 

The removal of the stellar contamination by interpolat¬ 
ing and subtracting the continuum emission worked well 
for all objects up to a certain brightness limit. Prob¬ 
lems occurred for saturated objects and their spikes as 
well as for those objects with diffraction rings. At the 
positions of these sources, the pixel values are not reli¬ 
able anymore. Therefore, the interpolation and thus the 
subtraction might partly or totally fail for these objects. 

For small areas with just a few unreliable pixel values 
we were able to clean the final reddening map by per¬ 
forming a spline interpolation. This effect appears for 
objects close to the saturation limit meaning that only 
a few pixels are saturated or the object is saturated just 


in one filter. 

For cosmic rays, dead and hot pixels, being these cases 
where just a single pixel is defective, we were able to 
clean the final map by replacing the affected pixels by 
the median value of their surrounding. 

For the very bright objects as well as for the whole clus¬ 
ter center, the area of saturated pixels was too large for 
a meaningful interpolation. Here we replaced the entire 
affected area by the median reddening calculated over 
an annulus around that area. For the map in Fig. [8] we 
additionally introduced a random Poisson noise to the in¬ 
terpolated pixel values (see Fig. III). Those values should 
thus be handled with caution. 

5. THE STELLAR REDDENING TOWARDS THE 
DIRECTION OF WD2 

In order to get a better insight into the physical pa¬ 
rameters of Wd2 we need to derive the variability of the 
reddening for the stellar population of the cluster. 

5.1. Transformation of the gas excess map to the 
stellar excess map 

After creating the pixel-to-pixel gas excess map E{B — 
V)g there was the need to translate this color excess into 
a stellar color excess E{B — V)* in order t o pro perly 
deredden the s tellar photometry. iCalzettil (jl997[ l and 
iCalzetti et al.l (120001) found a linear empirical relation 
between the stellar and the gas reddening {E{B — V)* ~ 
/ • E{B — V)g). In star forming regions, the dust in front 
of the stars is removed or destroyed by stellar feedback. 
Therefore, the stellar reddening is smaller than the red¬ 
dening of the gas. The interstellar medium (ISM) con¬ 
sists of a different composition (grain size, density, etc.) 
throughout the Galaxy. Therefore, this linear, empirical 
relation is only valid for a small range in color excess 
E{B — V) and needs to be determined for each cluster 
independently. 

We used the results of the ground based opti¬ 
cal spectral types of massive stars in Wd2 from 
Vargas Alvarez et al.l ()2013l) and lR,auw etlin ()2007l 
201 If) to obtain the bolometric m agnitude of each of 
those stars using the Tab 1.4 of iSoarke & Gallagli^ 
(|2007f) . Combining this information with the zero-age- 
main-sequence (ZAMS)E!3 derived from the isochrones 
of the PAdova and Triest e Stellar Evolut ion 
CodeF^ (hereafter: PARSEC, iBressan et al.l l2012r) we 
could obtain the theoretical intrinsic colors in the proper 
filter set of the WFPC2. We used this information to 
determine the color excess of the star combining the in¬ 
trinsic color and the color derived with the photomet¬ 
ric catalog of iVargas Alvarez et al.l (|2013l) {E[B — V)* = 
[B — V)* — {B — V)th) since our photometric catalog does 
not cover the H-band. We then used the proper transfor¬ 
mation of the photometry from the WFPC2 photometric 
system to the Johnson-Cousins system: 

A(Al - \2)g _ 

E{B - V), 

The ZAMS was derived, combining PARSEC isochrones with 
the information given on the webpage, at which age a star for a 
certain mass interval reaches the main sequence. 

http://stev.oapd.inaf.it/cmd 
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TABLE 6 

The transformation parameters 


Filter 

Filter system 

Xp [nm] 

a 

b 

F336iy 

WFPC2 

334.44 

0.87338 

2.37152 

F439VE 

WFPC2 

431.13 

0.99552 

1.12082 

B 

— 

445® 

1.00213 

0.94137 

V 

— 

551® 

1 

0 

F555W 

ACS 

536.10 

1.00699 

0.06883 

F555VE 

WFPC2 

543.90 

1.00311 

0.02704 

F814VE 

WFPC2 

801.22 

0.78049 

-0.57273 

F814VE 

ACS 

805.70 

0.77652 

-0.58020 

F125W 

WFC3IR 

1248.6 

0.40149 

-0.36821 

F160W 

WFC3IR 

1536.9 

0.28735 

-0.26382 


Note. — The parameters a and b for transforming the 
color excess of the stars E{Xi — A 2 ) to an excess E{B — V). 
The parameters a and b a re calculated with Eq. 3a and 3b 
from ICardelli et al.l 1119891 1. 

^The pivot wavelength of the B and V band are taken from 
IBinnev fc Merrifieldl II1998I . p.53, Table 2.1) 


Rv Ry 


as 


i>B 

Rv 


ay 


by 

Rv 


( 2 ) 


The para meters a and b are ca lculated with Eq. 3a 
and 3b from iCardelli et al.l (I1989D and given in Tab. [51 
Al and A 2 are the pivot wavelengths of the F4391E and 
F5551T filters, respectively. For the total-to-selective ex¬ 
tinction, we used the same value as the color excess map 
{Rv = 3.1). We should note here that even though we 
infer an Ry here in order to determine E{B — E)* even 
that we analyze Ry in Subsect. 17.11 A variability of Ry 
between 3.1 and 4.0 leads to an change in E{B — E)* of 

0.6%. In Table [71 we list the values of the gas excess 
from our gas excess map, as well as the calculated stellar 
excess, (Column 8-11) for each star. 

In order to get the transformation between the stellar 
excess and the gas excess we used all stars with spec¬ 
troscopy that have individual reddening values (e.g.: are 
not located in the center) and are main-sequence stars 
(class V without special line features). A linear, error 
weighted regression revealed the following relation be¬ 
tween the gas and the stellar extinction: 


E{B - V)^ = 0.4314 • E{B - V)g + 0.7400. (3) 

It is satisfactory to use a total-to-selective extinction 
of Ry = 3.1 for creating the color-excess map of the gas 
(see Sect. 11 since a different Ry will be compensated by 
different regression parameters in Eq. H and, therefore, 
leaves the overall outcome for E{B—V)i, unaffected. This 
results in an intersect greater than zero. Eq. O now 
allows us to translate the E{B — V)g gas extinction map 
into a map of the stellar extinction E{B — E)*. The 
overall median of E{B — V)g = 1.87 mag transforms to 
E{B — E)* = 1.55 mag. In Fig. [Hone can see the stars 
used to calculate the transformation. 

In Fig in the histogram of the distribution of the color 
excess E{B — V) is plotted for each of the 14199 objects 
in the catalog spatially covered by the color excess map. 
The gas excess distribution is shown in red while the 
stellar excess distribution is shown in blue. The asterisks 
mark the color excesses of the thirteen stars with spectral 
types, which we used to obtain the linear regression for 


transforming the gas excess into a stellar excess. It can 
be seen that the color excess of the 13 stars is equally 
distributed over the whole range in color excess of all 
objects in the catalog. Therefore, the linear approach 
is valid for the complete catalog. However, we caution 
that this transformation cannot be used to extrapolate 
to E{B — V)g = 0 mag or to very high E{B — V)g values. 



Fig. 14.— The 13 stars with spectroscopic data used to calcu¬ 
late the transformation between the pixel-to-pixel gas excess map 
E{B — V)g to the stellar excess E{B — V)* are shown here. The 
green line is the best fit linear regression including error weighting. 



E(B-V) 


Fig. 15.— The distribution of the color excess E(B — V) for each 
of the 14199 objects in the catalog, covered by the color excess 
map. The gas excess distribution is shown in red while the stellar 
excess distribution is shown in blue. The asterisks mark the color 
excesses of the thirteen stars with spectral types used to obtain 
the linear regression for transforming the gas excess into a stellar 
excess. 


5.2. The stellar color excess derived with UBV 
photometry 

An independent way of deriving the stellar redden¬ 
ing towards Wd2 is to use UBV two-color diagrams 
(TCD). Since our own photometric catalog lacks the 
H-band photometry, we used the photometric cata¬ 
log of iVargas Alvarez et al.l (I2013D . includ i ng th e spec¬ 
troscopy derived by iRauw et al.l (|2007L 120111) . and 
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TABLE 7 

Comparison to spectroscopic data 


RA 

J2000 

Dec 

J2000 

Spec. Type 

F555VE 

{F439W 

int,th 

- F555W) 
V.-A. 

E{B - V). 
TCD map 

* 

O ' map 

E(B 

map 

-V)g 

O ' map 

An. 

10:23:55.176 

-57:45:26.89 

04 V 

15.610 

-0.281 

1.776 

1.841 

1.526 

0.013 

1.822 

0.024 

b 

10:23:56:160 

-57:45:29.99 

04 V-III((f)) 

14.514 

-0.281 

1.451 

— 

1.488 

0.006 

1.734 

0.020 

b 

10:23:59.201 

-57:45:40.53 

07.5 V 

16.015 

-0.270 

1.666 

1.734 

1.629 

0.099 

2.061 

0.038 

b 

10:24:01.920 

-57:45:32.62 

08.5 V 

16.100 

-0.263 

1.489 

1.574 

1.686 

0.014 

2.193 

0.027 

ajb 

10:24:00.353 

-57.45:42.71 

08 V 

15.863 

-0.267 

1.615 

1.687 

1.528 

0.018 

1.827 

0.036 

b 

10:23:00.480 

-57:45:24.01 

04 V 

14.546 

-0.281 

1.485 

1.581 

1.514 

0.007 

1.794 

0.015 

b 

10:24:00.499 

-57.44:44.53 

B1 V 

15.587 

-0.254 

1.358 

1.441 

1.437 

0.018 

1.615 

0.030 

b 

10:24:00.713 

-57.45:25.42 

08 V 

15.466 

-0.267 

1.410 

1.503 

1.514 

0.015 

1.794 

0.023 

b 

10:24:00.816 

-57.45:25.87 

06.5 V 

14.866 

-0.276 

1.454 

1.549 

1.514 

0.012 

1.794 

0.015 

b 

10:24:00.979 

-57.45:05.50 

B1 V 

16.111 

-0.254 

1.336 

1.421 

1.473 

0.020 

1.699 

0.041 

ajb 

10:24:01.070 

-57.45:45.73 

09.5 V 

16.239 

-0.259 

1.482 

1.561 

1.576 

0.017 

1.938 

0.027 

b 

10:24:01.200 

-57:45:31.07 

03 V 

13.487 

-0.283 

1.513 

1.608 

1.543 

0.017 

1.862 

0.006 


10:24:01.392 

-57:45:29.66 

04 V 

14.084 

-0.281 

1.411 

1.515 

1.543 

0.006 

1.862 

0.007 


10:24:01.454 

-57.45:31.33 

03 V 

15.042 

-0.283 

1.403 

1.510 

1.543 

0.007 

1.862 

0.017 


10:24:01.524 

-57.45:57.06 

06 V 

15.053 

-0.278 

1.810 

1.868 

1.682 

0.079 

2.184 

0.017 

b 

10:24:01.610 

-57.45:27.89 

05.5 V 

14.921 

-0.278 

1.378 

1.482 

1.543 

0.008 

1.862 

0.027 


10:24:01.889 

-57.45:40.05 

09.5 V 

16.574 

-0.259 

1.656 

1.717 

1.543 

0.023 

1.862 

0.042 

a 

10:24:01.889 

-57.45:28.00 

08 V 

15.524 

-0.267 

1.510 

1.593 

1.543 

0.015 

1.862 

0.022 


10:24:02.064 

-57:45:28.01 

06 III 

14.453 

-0.278 

1.664 

— 

1.543 

0.010 

1.862 

0.009 


10:24:02.186 

-57.45:31.32 

09.5 V 

16.609 

-0.259 

1.601 

1.668 

1.543 

0.020 

1.862 

0.040 

a 

10:24:02.256 

-57.45:35.12 

04-5 V 

14.736 

-0.281 

1.498 

1.592 

1.543 

0.007 

1.862 

0.018 


10:24:02.304 

-57:45:35.53 

04.5 V 

13.893 

-0.281 

1.475 

1.571 

1.543 

0.006 

1.862 

0.008 


10:24:02.376 

-57:45:30.64 

03 V((f)) 

13.878 

-0.283 

1.618 

— 

1.543 

0.018 

1.862 

0.004 


10:24:02.414 

-57.45:47.11 

09.5 V 

16.637 

-0.259 

1.704 

1.760 

1.592 

0.112 

1.976 

0.048 

a,b 

10:24:02.448 

-57.44:36.13 

05 V-III 

13.155 

-0.280 

1.359 

— 

1.479 

0.052 

1.712 

0.003 

b 

10:24:02.518 

-57.45:31.47 

08.5 V 

15.790 

-0.263 

1.598 

1.671 

1.543 

0.014 

1.862 

0.033 


10:24:02.555 

-57.45:30.52 

08.5 V 

16.243 

-0.263 

1.649 

1.717 

1.543 

0.017 

1.862 

0.030 


10:24:02.604 

-57.45:32.26 

06-7 V 

16.118 

-0.276 

1.579 

1.661 

1.543 

0.016 

1.862 

0.039 


10:24:02.664 

-57.45:34.38 

03-4 V 

14.554 

-0.282 

1.528 

1.620 

1.543 

0.005 

1.862 

0.016 


10:24:02.789 

-57.45:30.05 

08 V 

15.963 

-0.267 

1.593 

1.667 

1.543 

0.017 

1.862 

0.035 


10:24:03.787 

-57.44:39.87 

09.5 V 

15.349 

-0.259 

1.380 

1.470 

1.543 

0.012 

1.862 

0.022 

b 

10:24:04.901 

-57.45:28.43 

04-5 V 

14.528 

-0.281 

1.606 

1.688 

1.511 

0.079 

1.787 

0.012 

b 


Note. — All stars with a determined spectroscopic type from IVargas Alvarez et al.l ||201.‘1H and IR.auw et ^ ] li20(y7Ll2Qnh . All 
WFPC2 photometry was adopted from fVargas Alvarez et al.l ll2013D . Columns 1—2 list the J2000 coordinates followed by the spectral 
type. Column 4 gives the WFPC2 F555W magnitude . Columns 5 gives the WFPC2 F439VF — F’SSSVF theoretical intrinsic colors, 
derived from (p.23, Tab.l.4. |Sparke fc Gallagheill2007l) in combination with the PARSEC ZAMS. Column 6 lists the stellar WFPC2 
^4391^ — F'bb'bW colors. Column 7 gives the stellar excess derived from the UBV TCD. Columns 8—11 list the stellar and gas 
excess E{B — V) and their respective uncertainties derived from our color excess map . Column 12 gives additional information: 
(a) stars have a proper ElbOVE and FSIAW photometry in our data and are shown in our CMD (Fig. |22l l . (b) individual reddening 
from the gas extinction map is available. 


IVargas Alvarez et alJ (j2013D . We used all stars with a 
spectral classification of luminosity class V, meaning that 
they are located on the main-sequence, and have no ab¬ 
normal line features (see Tab. 

Using the PARSEC ZAMS with the same procedure 
as described in Subsect. oi we were able to obtain the 
locus of the spectroscopically determined stars on the 
main-sequence in the UBV TCD. By shifting the indi¬ 
vidual stars in the TCD to their theoretical loci on the 
main sequence we could derive the stellar color excess 
A(F439W-F555W)* and A(E336W - F439W)*. Col¬ 
umn 7 of Tab. [7] lists the derived values transformed 
to the Johnson -Cousins photometric system with the ex¬ 
tinction law of iCardelli et al.l (1198911 (see Eq. 

The slope of the two-color excess in the TCD gives a 
ratio of E{U - B)/E{B -V) = 0.85 ± 0.033. Therefore, 
it is larger than the mean throughout the MW, whic h 
amountsto E{U—B)/E{B—V) = 0.72 (|Eitzgeraldlll970ll . 
This is consistent with the higher total-to-selective ex¬ 
tinction derived in Subsect. 17.11 The mean of the color 
excess of the stars derived via the TCD is E{B — U)* = 
1.62 ± 0.113. The median of the color excess of the stars 
derived with the color excess map was E{B — V)i, = 1.56. 
Comparing the individual values of E{B — U)* derived 


with the TCD with those of the color excess map reveals 
a mean difference of AE{B — U)* = 0.087 ± 0.129. We 
can say, within the errors, that our method of deriving 
the pixel-to-pixel gas excess map and transferring it to 
an excess map of the stars was successful and we can 
confidently use this map to deredden our photometry of 
the cluster members. 

6. THE COLOR-MAGNITUDE DIAGRAM OF WD2 
6.1. The morphology of the CMD 

In order to understand the stellar populations in our 
dataset we plotted the 7^81411^ — F160IF vs. F814IF 
color-magnitude diagram (GMD) of our catalog (see Fig. 
fTOl) . This CMD contains 7697 objects in total. Look¬ 
ing at the CMD it immediately becomes clear that we 
observe a composite population. There are two parallel 
sequences: a blue one at F8I4IF — FI60IF ~ 2-4 mag, 
and a redder one, at F814LF — F160IF ~ 4-8 mag (with 
approximate boundaries represented by the green lines 
in Fig. [m. The bluer population seems to consist of 
two branches, probably due to stars at different distances 
along the line of sight through the Galactic disk. The 
brightest objects can be seen at 7^81411^ ~ 14 mag. 
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The bluer sequence of the two clearly separated popu¬ 
lations can be explained with the foreground population 
of the Galactic disk towards Wd2 (see Sect. [ 3 ). This 
population is located closer to us and hence less red¬ 
dened. This makes those stars appear bluer in the CMP 
(see t he comparison to the Besangon model (IRobin et al.l 
1200^ in Fig. [H]). The red sequence of the CMD is most 
likely the cluster population of Wd2 (see Sect. 16.21) . 

Our data suffer from saturation of the very bright ob¬ 
jects (see Sect. |31), and for this reason we do ’’lose” 
stars in the CMD at the bright end. Due to the de¬ 
tection method in the WFC3-IR filters (see Sect. [2]) 
saturation occurs only in the optical filters. Compar¬ 
ing our objects with detections in at least two filters 
with more shallow photometric studies from the liter¬ 
ature ([Vargas Alvarez et al.l 120131) we conclude that we 
miss 16 luminous stars, at the bright end of our photom¬ 
etry. 



Fig. 16.— The FS14W vs. (_F814VF — F160W) CMD of the Wd2 
region. The plot contains 7 697 objects. On the right representative 
photometric error bars are shown in bins of 1 mag. The green 
lines represent the separation between the two parallel sequences 
of likely foreground stars (blue sequence) and point sources likely 
associated with Wd2. 


6.2. Selection of likely RCW 49 member stars 

In order to study the properties of Wd2 and the spatial 
morphology of the stellar content we need to extract a 
sample of probable RCW 49 member stars. Therefore, 
we need to clean the catalog from foreground objects (see 
Sect. O. 

The RCW 49 member selection was performed using 
the F8141F - F160W vs. F8UW CMD of our catalog 
(see Fig. nil) which contains a total number of 7 697 ob¬ 
jects (see Sect. El). As a first step, we divided the data 
in likely foreground population stars (5 604) (marked 
by black points) and the likely RCW 49 member stars 
(2 093), indicated by red points (see Fig. [T71). Distin¬ 
guishing these two populations by defining a diagonal 
line was possible due to the clear separation of the two 


sequences and even at the faint end of the CMD the pho¬ 
tometric uncertainties are still small enough to clearly 
separate the sequence formed by unrelated Galactic disk 
objects and RCW 49 members. 

To test whether the black objects are indeed contam¬ 
inant stars of the Galactic plane, we simulated the MW 
disk in the direction of Wd2 with the Model of stellar 
popula tion synthesis o f the Galaxy, the Besangon 


model (IRobin et al.l 1200311. This model adoots the ex- 

tinction law of Mathis 

(1 19901). It lacks the spiral arm 

structure ([Robin et al.l 

200311 ■ but since we are just sam- 


pling one line of sight, this should not be a prob¬ 
lem. We simulated the color-magnitude locus of the 
expected stars in the Galactic plane up to the clus¬ 
ter distance of 4.16 kpc across a FOV corresponding in 
size to our survey area (~ 21 arcmin^). The simulated 
stellar distribution is visualized by the green points in 
Fig. [m We transformed the model colors and bands, 
which are by default in the Johnson- Cousins photo¬ 
metric system (iJohnson &: Morg^ll953ll . into the HST 
photometric system with the SYNPHOt/calcphot rou- 
tinj^ ljhaidler et ^1200511 for consistency with the data. 
It looks like the green points are centered between the 
two slightly separated populations of the foreground, 
most likely originating from a non-uniformly distributed 
ISM in the line of sight towards Wd2. The green points 
do indeed nicely overlap with the blue diagonal sequence 
of suspected foreground stars, supporting our assump¬ 
tion. Since we do not know the distribution of the fore¬ 
ground reddening we cannot reproduce the exactly ob¬ 
served distribution of the foreground populations. 



0 2 4 6 8 10 

F814W-F1 60W 


Synphot is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA. 
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Fig. 17.—The F814ty vs. (F814VF-F160VF) CMD of the 7697 
point sources from the Wd2 photometric catalog. The red points 
represent likely members of the RCW 49 region (including cluster 
members). The black sources are probable foreground stars. The 
green points indicate the location of the synthetic objects from 
the Besangon model in the direction of Wd2 up to the cluster 
distance of 4.16 kpc falling within a FOV corresponding in size 
to our survey area (~ 21 arcmin^). Their superposition with the 
black points suggests that the latter are less reddened foreground 
stars. The black solid line marks our arbitrary separation between 
cluster and field stars. 


7. THE PHYSICAL PARAMETERS OF WD2 

Using our high resolution dataset in combination with 
the 2D color excess map we are able to determine the 
physical properties of Wd2. 

7.1. The total-to-selective extinction towards the 
direction of Wd2 

In order to correct for the high differential reddening, 
visible in the map of the color excess E{B — V)g (see 
Fig. m), we need to get a better insight in the total-to- 
selective extinction ratio Rv'- 


Rv 


Ay 

E{B - V) ’ 


(4) 


where Ay represents the total extinction in the visible. 

Since the distance modulus for Wd2 is not yet fixed, we 
used the TCD to derive a value for Ry, since TCDs have 
the advantage of being distance independent. To derive 
Ry, we had to find the best fit of the ZAMS to the likely 
main-sequence stars of our photometric catalog. There¬ 
fore, we used all stars brighter than F814IF > 18 mag 
of the RCW 49 members (red points in Fig. [18] and Fig. 
[TOl selected through the CMD of FigfTTl). We used three 
different TCDs: F5551F-F8141F vs. E8UW-E125W, 
F555W - F8UW YS. F125W-F160W, and {FA39W- 
F555W)yYFPC2 vs. (F5551F — F814lW)wfpC2- Here 
again, we also used the WFPC2 data from the photomet¬ 
ric catalog of|Wrg^3I^siiIiE^ The PARSEC 
ZAMS were used for Solar metallicity. 

In order to deredden the stars onto the locus of the 
ZAMS in the TCD s we used the extinction law of 
iCardelli et al.l (|1989ll : 


^A= (^) ■A(H-U)*.Rv. (5) 

A\/Ay is also i?v-dependent: 


A 

Ay 


— flA + 


Ry 


( 6 ) 


The parameters a and b are wavelength dependent and 
given in Tab. [HI Combining Eq. [Sjand Eq. jSjwe can cal¬ 
culate the absolute extinction A\ in dependence of Ry. 
In order to get the best fit for the three TCDs we needed 
to apply values of Ry = 3.81 for F555IF — F8\AW 
vs. F814IF - F125VF (see Fig. |TH|), i?y = 3.96 
for F555IF - F814IF vs. F125IF - F1601F (see Fig. 
nH), and Ry = 4.08 for (F439IF - FbbbW)wFPC2 vs. 
(F5551F - F81AW)wfpC2 (see Fig. |20|). In the cor¬ 
responding Figures [181 [TH and |20| we show the stellar 


population including the ZAMS and the reddening vec¬ 
tor. The red points represent the likely main-sequence 
stars, the blue diamonds represent the objects for which 
spectroscopy is available. The discrepancies at the red 
end of the IR filters in TCD to the remaining objects 
(black points) originate from the very numerous PMS 
population appearing redder due to remaining circum- 
stellar dusty envelopes. Due to the high number of spec¬ 
troscopically observed stars in Fig. |20l we only used 
those objects for fitting. The ZAMS fitting gives us a 
mean value of Ry = 3.95 ± 0.135. Spectroscopic obser- 
vations of the most ma s sive s tars (|Vargas Alvarez et al.l 
120131: iRauw et al.]l2007L 1201111 and a photometric study 
of lHiir et al.l (|2015fl show that Wd2 has a total to selec¬ 
tive extinction in the visual {Ry = Ay/E{B — V)) of 
about 3.64-3.85, what agrees with our value within the 
uncertainties. 



Fig. 18. — The F555tF — F814VF vs. F814VF — F12bW 
TCD of the Wd2 photometric catalog. The green line represents 
the PARSEC ZAMS. We show in red the reddening vector for 
E{B — V)* = 1.55 mag and Ry = 3.81. The red points represent 
the likely main-sequence stars (all RCW 49 members brighter than 
F814W > 18 mag). The discrepancies at the red end originate 
from the rich PMS population appearing redder due to the per¬ 
sisting dusty envelopes around the stars. The blue diamonds mark 
the stars with spectral types from (IVargas Alvarez et al.ll2ni3lb 


7.2. The extinction correction of our cluster CMD 

The high resolution color excess map of the gas E{B — 
V)g, was transform ed t o the stellar color excess E{B — 
U)* (see Subsect. 15.1|) and used to correct the CMD 
for differential reddening (see Sec. 14.4p . In other words, 
we used the stellar color excess at the position of the 
stars to subtract the foreground reddening. For the color 
excess E{B — U)* at each stellar position, the median 
E{B — U)* within three sigma of the mean PSF was 
calculated and then translated into the total extinction 
A\ at the pivot waveleng th A of the HST filte rs, using 
again the extinction law of iCardelli et al.l (I1989D (see Eq. 
|5|). Ax/Ay is given as 0.62303 for F8141F and 0.21756 
for F1601F. 

In Fig. |2T] we show the CMD of the Wd2 members 
as observed (left) as well as the individually dereddened 
CMD (right). The red arrow shows the overall median 
reddening of E{B — U)* = 1.55 with Ry = 3.95 (see 
Subsect. 17.11 for determining Ry). We should note here 
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(FI 25W-F1 60W)o 


Fig. 19.— The F555VF - F814VF vs. F125VF - F160VF TCD 
of the Wd2 photometric catalog. The green line represents the 
PARSEC ZAMS. In red we show the reddening vector for E{B — 
F)* = 1.55 mag and Ry = 3.96, shortened by a factor of 4 for 
plotting reasons. The red points represent the likely main-sequence 
stars (all RCW 49 members brighter than F814VF > 18 mag). The 
discrepancy at the red end originate from the rich PMS population 
appearing redder due to the persisting dusty envelopes around the 
stars. The blue diamonds mark the stars with spectral types from 
(IVargas Alvarez et al.ir2013IV 



Fig. 20.— The (F439IT — F555W)wfpC2 vs. (F555VF — 
F8UW)iyFPC2 TCD of the WFPC2 photometric catalog of 
IVargas Alvarez et al.l II2013I '). The green line represents the PAR¬ 
SEC ZAMS. The reddening vector for F{B — V)* = 1.55 mag and 
Ry = 4.08 is shown in red. The blue diamonds mark the stars 
with spectral types from ll Vargas Alvarez et al.ll2013ll . 

that objects lying in regions where the color excess map 
needed to be in terpolated or substituted by the median 
(see Sect. 14.4L the uncertainty of the dereddening is 
bound to the uncertainty of the estimation of the values 
of the saturated pixels. Consequently, individual dered¬ 
dening means using the individual stellar excess for each 
star and the extrapolated/substituted value for the stars 
where the color-excess map has corrupt pixels. As one 
can see from the color representation of the density in 
the plot, a tightening of the CMD was achieved. Nev¬ 
ertheless, the spread cannot be reduced to less than the 
intrinsic age spread of ^ 2 Myr. Also effects like rotation 
and and binarity play a role. PMS stars have their own 
circumstellar effects, in addition to the general differen¬ 


tial reddening, causing broadening of the CMD. 

7.3. The distance and age 


The distance to Wd2 is still subject of debate in 

the literatur e. _ Pu blished values range from 2.8 kpc 

(|Rauw et al.l 120071 : lAscenso et al.l l20()7t iCarraro et al.l 


I2013D. to 4.16 kpc (IVargas Alvarez et al.ll2013l). 5.7 koc 

(IPiatti et al.lll998D. 6.4 k 

pc (Carraro & Munarill2004D to 

8 koc dRaiiw et al.ll2007l 

201 1|). Discrepancies also exist 


in the published ages: an upper age limit of 3 Myr was 
inferred for th e whol e clu ster and 2 Myr for th e core by 
lAscenso et al.l (I2007D and ICarraro et al.l (I2013D . 


In this paper, we will attempt to derive tighter con¬ 
straints on distance and age from our high resolution, 
high sensitivity multi-wavelength observation. 

In order to determine the distance and age of Wd2 
we used the method of over-plotting isochrones to our 
F8UW vs. {F8UW - F160VP) CMD. We especially fo¬ 
cused on the turn-on, since here the isochrones are not 
degenerate in age and most distinguishable between each 
other, compared to the main sequence. We used PAR¬ 
SEC isochrones for 0.5-2.0 Myr and Solar metallicity. 

In order to determine the best distance, we used the 
already dereddened CMD (see Fig. [521 rig ht panel). We 
corre cted the chosen PARSEC isochrones (iBressan et al.l 
l2012t) for the distance modulus. We applied different 
values for the distance and assessed the best fit to both 
the locus of the turn-on region and the PMS population. 
This corresponded to a distance of d = 4 kpc and is in 
good agreement with the spectroscopic observations of 
IVargas Alvarez et al.l (I2013D . Their distance is based on 
the physical constraints derived from the spectral types 
of the most massive stars in Wd2. The method they 
adopted is more accurate than our empirical fit, and, 
therefore, we will assume for all the subsequent calcula¬ 
tions their value of d = 4.16 kpc. 

Taking into account the magnitude range of the turn¬ 
on we can conclude that the age for the whole cluster is 
in the range of 0.5-2.0 Myr. We provide a more detailed 
and spatially differentiated analysis of the age in Sect. 
[51 The uncertainty in the total-to-selective extinction 
ARy ± 0.135 results in an age uncertainty of ~ 0.5 Myr 
(see dashed-dotted lines in Fig. [52]) assuming that the 
luminosity of the stars is constant. The impact on the 
distance is ±0.01 kpc and, therefore, negligible. 

8. THE TWO CLUMPS OF WESTERLUND 2 

8.1. The spatial distribution of the stellar content 

As can be seen from Fig. |21 our observations cover 
an area ^ 4 arcmin wide, which, assuming a distance of 
4.16 kpc, corresponds to a linear scale of 4.8 pc. This re¬ 
gion covers the Wd2 cluster and the surrounding clouds, 
and enables us to analyze the spatial distribution of the 
stars. Due to the very complex cloud structure and the 
resulting variable reddening (see Sect. H as well as the 
different degrees of crowding, we expect that the com¬ 
pleteness of our dataset differs as a function of posi¬ 
tion. In order to avoid the bias due to incompleteness ef¬ 
fects we applied luminosity cuts based on the histogram 
of the magnitude distribution (see Fig. ini)- In other 
words, we limited our analysis to the range where the 
catalog is mostly complete. We used all objects with 
F8141U < 22.75 mag and F160IU < 18.00 mag. This 
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(F81 4W-F1 60W) (F81 4W-F1 60W)o 


Fig. 21.— Left: The observed ^814144 vs. (F814VF — FlbOlF) CMD of the likely cluster members. Right: The F814VF vs. {F^lAW — 
FIQOW) CMD of the of the individually dereddened cluster stars. The red arrow marks the reddening vector corresponding to the median 
E{B — y)* = 1.55 mag with Ry = 3.95. The colorbar represents the number of objects per bin (binsize: 0.1.2 mag). The range in color is 
the same in both plots. 



(F814W-F160W)o 

Fig. 22.— The extinction corrected F814VF vs. (F814VF — 
FIQOW) CMD of the 5 418 point sources left the Wd2 member se¬ 
lection. The colored lines represent the PARSEC PMS isochrones 
for 0.5—2.0 Myrs best fitting the turn-on of the main sequence. 
The error bars represent the typical photometric uncertainties in 
one magnitude bins. We applied an Ry = 3.95 with a distance of 
4.16 kpc. The red diamonds mark the stars for which spectroscopic 
data are available (see Tab. [7|l. The dashed-dotted lines represent 
the maximum possible shifts of the isochrones for ARy ± 0.135 
assuming a constant luminosity of the stars. 

leaves us with 2 610 point sources. A detailed assessment 
of the incompleteness through artificial star experiments 
will be performed and published in a subsequent paper. 

Using this reduced sample we created a 2D plot of the 
stellar density, taking into account all objects within a ra¬ 
dius of r < 10.8" around each star. This area was chosen 
such that there were always enough objects within the 
area and the resolution was high enough to resolve the 
spatial structure of the stellar content. The result of this 
approach can be seen in Fig. [531 Each dot represents one 
star. Fifteen linear spaced lines representing isodensity 
contours are over-plotted. The color indicates the num¬ 
ber density specified by the color bar [number/arcmin^]. 

In the resulting surface density plot the cluster appears 
divided into two separate clumps; the main cluster of 
Wd2 and a smaller clump to the north (northern clump). 
This structure was also pointed out and named northern 
clump bv iHur et al.l (|2015D . iVargas Alvarez et al.l (|2013[1 


also noticed a ” a secondary concentration of stars 45" to 
the north”. 

In order to investigate the properties of and possible 
differences between the two clumps we defined a circular 
area around each clump containing stars that we consid¬ 
ered to belong to each clump. We removed the uniform 
distribution of stars, visible outside the cluster, by deter¬ 
mining a threshold density. This density was determined 
in a control field well outside the central cluster. All 
objects 3cr above this stellar density were included. 

We then plotted histograms in right ascension and dec¬ 
lination of all remaining objects (Fig. (Ml) . Fitting Gaus- 
sians to the distribution and determining the FWHM 
allowed us to estimate the extent of the area where the 
two clumps dominate. The fits were always made under 
the condition that the peaks in right ascension and dec¬ 
lination for the same clump needed to be similarly high 
in order to be consistent. 

For the declination direction we used a combination of 
two Gaussians (dashed-dotted-line) (see left of Fig i24)) . 
The red curve fits the maximum of the histogram dis¬ 
tribution with a center coordinate of —57°45'31"7 and a 
FWHM of 26.83". This coordinate marks the locus of 
the main body of Wd2. The blue curve fits the north¬ 
ern clump with a center coordinate of —57°44'34"2 and 
a FWHM of 19.17". 

For the right ascension direction, we needed to subdi¬ 
vide the data since the two clumps are located at very 
similar right ascensions. Therefore, we divided the cat¬ 
alog at Dec=—57°44'51"0. We plotted the histograms 
for each part (see right side of Fig. [M]). Stars associ¬ 
ated with the main body of Wd2 are plotted in red and 
those with the northern clump in blue. The red curve 
fits the main body of Wd2 with a center coordinate of 
RA= 10^24’"01".63 and a FWHM of 24.76". The blue 
curve fits the northern clump with a center coordinate of 
RA= 10'‘24’"01®.63 and a FWHM of 24.76". The Gaus¬ 
sian fit of the northern clump is higher than the actual 
distribution due to the fact that for both declination and 
right ascension the peaks needed to be of similar height. 
Since we do not have any completeness information, yet, 
we assumed a circular shape for the clumps. We will be 
able to determine more accurately the spatial structure 
of those clumps after having assessed the completeness 
of our photometry. This will be included in a subsequent 
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Fig. 23.— The surface number density of the Wd2 stellar content. Each dot represents one single star. The color bar indicates the density 
in number arcmin“^. The contours show the density in 15 linear steps. The outer boundaries of our survey are indicated by thick black 
lines. The approximate radii of Wd2 = 24.10^^) and the northern clump {vnc = 18.96^^) determined by the FWHM of a Gaussian 

distribution are shown by two large circles. 


paper. 

We calculated the radii of the clumps as the mean of 
the radii along the declination and right ascension direc- 
tions, and found a radius of rwd 2 = 24.10" 0.5 pc) for 

the stars associated with the main body of Wd2 and a 
radius of r„c = 18.96" (~ 0.4 pc) for the northern clump. 
These radii are plotted as circles in Figl23l These radii 
merely serve to select stars that mainly belong to one of 
the two clumps. A thorough analysis of the density pro¬ 
files (including incompleteness corrections) are deferred 
to a subsequent paper. 

8.2. The properties of the two clumps 
8.2.1. The morphology 

In Fig. [5n]one can see a gray-scale image of the central 
region containing the two clumps in the FSIAW band 
(left panel) and the reddening map of the same area 


(right panel). The stellar density contour lines from Fig. 
ESI and the circle marking the central regions of the two 
clum ps are over-plotted. Our spatial analysis (see Sect. 
18.1|) shows that the surface density of the stellar content 
is higher in the clumps than in the surrounding area. The 
reddening map shows an extended dust filament superim¬ 
posed on the northern clump. The de-reddened CMDs of 
both clumps (see Fig. 126)) are not shifted in color which 
leads to the conclusion that this dust filament is located 
in front of the cluster. In the main body of Wd2 the con¬ 
tour lines follow the gas excess towards the North-West, 
so it looks like the stars there already cleared some of the 
gas. Towards the South and South-East the gas excess 
increases and the number density of the stellar content 
decreases quite rapidly. 

8.2.2. The luminous stellar population 
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-57:48:00 -57:46:48 -57:45:36 -57:44:24 -57:43:12 10:24:31,2 10:24:19.2 10:24:07.2 10:23:55.2 10:23:43.2 

Dec (J2000) RA (J2000) 

Fig. 24.— Histograms of the number of stars as a function of spatial position (expressed in right ascension and declination) with a binsize 
of 5 arcsec. Left: Histogram in the declination direction with a superimposed curve fitting two Gaussians (dashed-dotted line). The red 
curve is a fit of the maximum of the histogram distribution and marks the center of the main concentration of Wd2 which we find to be 
at a declination of —57°45'3l''7 and to have a FWHM of 26.83^^ The blue curve fits the northern clump with a center coordinate of 
-57A44'"34S2 and a FWHM of 19.17". Right: Histograms along the right ascension direction for the main body of Wd2 (red) and the 
northern clump (blue). The red curve fits the main body of Wd2 with a center coordinate of 10:24:01.63 and a FWHM of 24.76". The 
blue curve fits the northern clump with a center coordinate of 10:24:01.63 and a FWHM of 24.76". 


Both clumps of Wd2 coutain multiple OB stars. Most 
of them are located iu the maiu body of Wd2. All OB 
stars with available spectroscopic data are saturated iu 
our data. 

There are 3 stars with spectroscopic data available 
fromiRauw et al.l (I2007LI2MTI) . aud iVargas Alvarez et al.l 
(|2013[) within the northern clump and 19 within the main 
body Wd2. The northern clump contains a B1 V, 05 V- 
III, and an 09.5 V star. The 19 stars in the main body 
of Wd2 are all 0-type stars with a spectral type ranging 
from 09.5 V to 03-4 V. 

8.2.3. The color magnitude diagrams 

In the next step we take a look at the dereddened 
CMDs (see Fig ESl) of the main body of Wd2 and the 
northern clump. We are assuming the two clumps to be 
at the same distance. We again use the whole PMS pop¬ 
ulation, taking into account all stars in the respective 
areas. This yields 969 objects for the main body of Wd2 
(left panel of Fig. 1^51) and 474 for the northern clump 
(right panel of Fig. [5^. PARSEC isochrones for ages of 
0.5-2.0 Myr are overplotted. 

As already seen in the images (see Fig. [25]) and in 
the spectra the main body of Wd2 contains by far more 
luminous (> 14 mag) stars than the northern clump, 
both along the clearly visible turn-on of the PMS, as 
well as on the main sequence. In the northern clump 
the turn-on is barely visible and there are almost no 
stars on the main sequence up to our saturation limit of 
^ F814IF = 13 mag. One reason of this apparent lack of 
more massive stars may well be an effect of the stoch as- 
tic sampling of the mass function (e.g., IBaked 119^ of 
the relatively sparsely populated northern clump. For a 
better comparison of the CMDs of the two clumps we re¬ 
duced the number of objects in the main body of Wd2 to 
the same number as in the northern clump by a random¬ 
ized removal of objects. In Fig. ETlone can see the such 
reduced CMD of the main body of Wd2 in black and the 
CMD of the northern clump in red. Besides the fact that 
the turn-on is more apparent in the main body of Wd2 
there is no visible difference. The scatter of stars located 
on the turn-on is similar for the main body of Wd2 and 


the northern clump. This leads us to suggest that both 
clumps formed at the same time, but what we observe as 
the main body of Wd2 presumably formed from an ini¬ 
tially more massive and possibly denser gas clump. The 
mass function (which is left for a subsequent paper) will 
reveal whether there is a substantial difference between 
the two clumps and their formation process. 

We note that the lack of 03/04 stars on the main se¬ 
quence of the northern clump may imply a slightly higher 
age, but may also just be a consequence of stochastic 
sampling effects as pointed out before. In any case, for 
such young ages isochrones do not permit us to constrain 
the age more accurately, and within the accuracy of our 
age estimates we may consider the two clumps to be co¬ 
eval. 

Our CMDs show a luminosity spread among the PMS 
stars in the turn-on region that, if taken at face value, is 
compatible with an age spread of 2 Myr. However , 
considering the findin gs of, e.g.. iBaraffe et ap (1200911 ; 
iHosokawa et al.l (1201 III , and iSoderblom et al.l ^ 20141) on 
the impact of accretion and differences in the survival of 
circumstellar disks, we conclude that the adoption of ~ 1 
Myr as cluster age is reasonable. This is also supported 
by the fact that the locus of the turn-on region is very 
sensitive to t he age (^ 0.33 mag M yr~^ at 3 Myr), which 
is shown by iCignoni et al.l (|2010D for the young cluster 
NGC 346. 


9. SUMMARY AND CONCLUSIONS 


In this paper we described the reduction of our re¬ 
cently obtained HST (ACS and WFC3) data of the very 
young star cluster Wd2. We created a photometric cata¬ 
log of 17121 objects in six different filters, four wide-band 
(F555IU, F814IU, F125W, F160IU) and two narrow- 
band filters (F658IV, F128iV). 

Using the two narrow-band filters sampling the Ha and 
Pa/3 line, we created the first high resolution pixel-to- 
pixel map of the color excess E{B — V)g of the gas in 
Wd2. This map has a median color excess E{B — V)g = 
1.87 mag and reveals considerable differential reddening 
across the region. A radial analysis of the color excess for 
the four quadrants (NW, NE, SW, SE) shows an almost 
constant increase from ~ 1.86 mag to values of about 
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Fig. 25.— Left: Gray-scale image of the main body of Wd2 and the northern clump in the FSIAW band. Right: Reddening map 
shown on the same spatial scale and the same intensity scale as Fig fSl B oth panels show the same central areas. The stellar density 
contours are over-plotted as black (left) and white (right) lines (see Fig. I23II . The circles denote the areas dominated by stars belonging to 
the main body of Wd2 and to the northern clump, respectively. 


2.15 mag to the Southwest and ~ 2.0 mag to the South¬ 
east. To the North we have a decrease of the excess down 
to ~ 1.7 mag due to missing gas and dust. The median 
error due to the Poisson statistics of the electron count on 
the detector amounts to 0.032 mag. We estimate a con¬ 
tamination of the continuum in f 1251V origin a ting f rom 
the Pa/3 line of ~ 22%. Following iPang et al.l (120111) we 
infer a contamination of 10-15% in the broad band filters 
due to other emission lines, which leads to an uncertainty 
of 0.1-0.15 mag on E{B — V)g. 

Using the spectroscopy catalog obtained by 
Vargas Alvarez et al.l (1201311 : IRauw et'H] (1200711 . and 
Rauw et al.l (1201 ill we translate the pixel-to-pixel gas ex¬ 
cess map E{B — V)g into a stellar excess map E{B — V)i,. 
The median of this map is E{B — V)i, = 1.55 mag which 
is consistent with the mean value for the stellar excess 
derived with the UBV TCD using the photometric 
catalog of IVargas Alvarez et al.l (1201311. __ 

Taking the PARSEC ZAMS (jBressan et al.ll2012ll and 
fitting them to three different TCDs (F5551V — F814tU 
vs. F8141V - E12bW (see Fig. [HI), F5551V - F814tU 
vs. U125TV - EimW (see Fig. dll), and (F4391V - 
U5551V)wi^pC2 vs. (F 5551 V-F 8141 V)wfpC 2 (see Fig. 
dH)) allowed us to derive a value for the total-to se¬ 
lective extinction of Rv = 3.95 ± 0.135. This is in 
good agreement with the studies of iVargas Alvarez et al.l 


(IMl : iRa^ et al.l (l2?i(T7l[20TTll and IHur et all (IM^ . 
who showed that Wd2 has a total to selective extinc¬ 
tion in the visual {Ry = Ay/E{B — V)) of about 3.64- 
3.85. The value of Ry lies well above th e mean of the 
inters tellar medium (ISM) of Ry = 3.1 (jCardelli et al.l 
1198911 . This indicates that the dust grains along the 
line-of-sigh t towards Wd2 are larger th a n normally in 
the MW (iMathis fc WallenhorstI Il981l : iMathisI 119901 : 
iCardelli et al.l I1989H . This fact was also indicated by 
iVargas Alvarez et al.l (|2013ll . This is consistent with the 
higher than normal s lope of the E(U — B)/E{B — V) 
TCD of 0.90 ± 0.035 (lTurneilll973L ll^ . 

In combination with the derived color excess map, the 
Besangon model, and the U814IU vs. (F814IU—F160IV) 
CMD we were able to separate the cluster stars from 
the foreground contamination in the region of the giant 
H II region RCW 49. The CMD shows two well sepa¬ 
rated populations, where the bluer one consists of fore¬ 
ground stars representing the field population towards 
the Carina-Sagittarius spiral arm, while the redder one 
can be considered to be part of RCW 49 cloud containing 
the Wd2 cluster. We dereddened the selected RCW 49 
population using the E{B — U)* map which was trans¬ 
lated into the total extinction at the pivot wavelength of 
the H ST filters using the extinction law of lCardelli et al.l 
(Il989f) . Using PARSEC isochrones we estimated a 
distance of d ~ 4 kpc, which is in good agreement 
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(F81 4W-n 60W)o (F814W-F160W)(, 

Fig. 26.— The _F814PF vs. {FSIAW — F160W) CMDs of Wd2 (left, 969 objects) and of the northern clump (right, 469 objects). The 
colored lines represent the PARSEC PMS isochrones for 0.5-2.0 Myr fitting the turn-on of the main sequence. The red diamonds in the 
left panel represent those objects for which spectroscopic data is available (see Tab. 0 



(F814W-F160W)o 


Fig. 27.— The _F814VF vs. (_F814VF—-F160VF) CMDs of the main 
body of Wd2 (black points) and the northern clump (red points). 
The number of objects in the main body of Wd2 was scaled to 
the number of objects in the northern clump (469) to facilitate 
the comparison. The colored lines represent the PARSEC PMS 
isochrones for 0.5-2.0 Myrs fitting the turn-on regions of the PMS 
stars towards the main sequence. 


with t he more accurate result of iVargas Alvarez et alJ 
lj2013l l oi d = 4.16 kpc, which we ultimately adopted. 
With this distance Wd2 i s located at the sanie dis- 
tance (|Gennaro et al.l 120111: iKoumoia fc BonanosI 1201211 


as Westerlund 1 (|Westerlundl 1196111 . the most massive 
open star cl uster known in our Gal axy, and is closer than 
NGG 3603 (jHaravama et al.lF2008ll . 

The spatial density map (see Fig. 123)) of the stellar 
content of Wd2 reveals the splitting of Wd2 into two 
clumps (denoted as the main body of Wd2 and the north¬ 
ern clump). We infer the approximate area in which ob¬ 
jects can be assumed to belong primarily to the respec¬ 
tive clumps. This gave us projected radii (FWHM of the 
Gaussians) of rwd 2 = 24.10" (~ 0.5 pc) and rnc = 18.96" 
(^ 0.4 pc) for Wd2 and the northern clump, respec¬ 
tively. This leads to a projected area of ^ 0.5 arcmin^ 
for the main body of Wd2 and ^ 0.3 arcmin^ for the 
northern clump. Including all stars with F814IF > 
22.75 mag and F160IF > 18.00 mag, the mean stel¬ 
lar density of the main body of Wd2 is approximately 
2 000 stars arcmin“^, while the northern clump has a 


density of roughly 1600 stars arcmin“^. The stellar num¬ 
ber density is therefore ~ 20% smaller in the northern 
clump. 

A comparison of the CMDs for the two clumps 
showed no age difference. PMS isochrones indicate 
for both an age between 0.5-2.0 Myr. This leads us 
to conclude that both clumps were born at roughly 
the same time (l-OlJ;® Myrs ago). Therefore, Wd2 
has the same age as other very young s tar clusters 
like N GC 3603 in th e Carina spira l arm (jPang et al.l 
l201.3fl . Trum pler 14 (iCa rraro et ap 1200411 in the Ca- 
rina Nebula (ISm ith fc B rooks! ~ 08ll. R 1 36 in the LMC 
( Walborn fc Blades! ll99'11 :~ Sabbi et al.l 1201211 . Arches 
( Figer et al.l l2002t iFiged 12005 1. and is younger than 
Westerluiid 1 ( Clark et al.l l2005t iGennaro et al.l 120111 : 
iLim et ani2013[ l. The main body of Wd2 is more massive 
than the northern clump and shows a well defined PMS 
turn-on region. The main body of Wd2 contains at least 


19 O - type stars (|Vargas Alvarez et al.ll2013t iRauw et al.l 
120071 [2f)Tlh while the northern clump just contains one 1 
B and 2 0-type stars leading to a significant stellar mass 
difference. 

Such subclusters appear to be a common occurrence 
in massive star Terming regions, comprising between 2 to 

20 objects (see lKuhn et al.l[20l4 and references therein). 
Wd2 thus shows a modest example of such a substruc¬ 
ture. 

In a forthcoming paper, we will carry out complete¬ 
ness tests and determine the present-day mass func¬ 
tion as well as the m ass of both clumps. iNota et ^ 
(IM5 . in preparationll describes the morphology of the 
region, the detection of a rich population of pre-main- 
sequence stars in the cluster, and a number of unusual 
protostellar objects observed in the Wd2 region. 


Based on observations made with the NASA/ESA 
Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. These observations are 
associated with program #13038. 
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